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Abstract 

The molecular mechanism of slip at the interface between polymer melts and weakly attractive 
smooth surfaces is investigated using molecular dynamics simulations. In agreement with our 



(N 



t4-H ' previous studies on slip flow of shear-thinning fluids, it is shown that the slip length passes through 

O ' 
"^ ■ a local minimum at low shear rates and then increases rapidly at higher shear rates. We found that 






at sufficiently high shear rates, the slip flow over atomically flat crystalline surfaces is anisotropic. 
It is demonstrated numerically that the friction coefficient at the liquid-solid interface (the ratio 



Q ' of viscosity and slip length) undergoes a transition from a constant value to the power-law decay 

o : 

as a function of the slip velocity. The characteristic velocity of the transition correlates well with 
►^ , the diffusion velocity of fluid monomers in the first fluid layer near the solid wall at equilibrium. 

^f-\ I We also show that in the linear regime, the friction coefficient is well described by a function of a 

in 

'!;;j- ' single variable, which is a product of the magnitude of surface-induced peak in the structure factor 

L»^ , and the contact density of the adjacent fluid layer. The universal relationship between the friction 

coefficient and induced fluid structure holds for a number of material parameters of the interface: 

fluid density, chain length, wall-fluid interaction energy, wall density, lattice type and orientation, 

sj . thermal or solid walls. 
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PACS numbers: 68.08.-p, 83.80.Sg, 83.50.Rp, 47.61.-k, 83.10.Rs 



I. INTRODUCTION 



The study of fluid transport through micro- and nanochannels is important for biotech- 
nological apphcations and energy conversion processes [l| . The precise control and manipu- 
lation of fluids in systems with large surface-to-volume ratios, however, require fundamental 
understanding of flow boundary conditions. Fluid velocity profiles can be significantly mod- 
ified in the presence of slip at a solid surface. The degree of slip is quantified by the slip 
length, which is defined as a distance between locations of the real interface and imaginary 
plane where the extrapolated tangential velocity component vanishes. It was shown by 
numerous experimental studies that the main factors affectingslippage at the liquid-solid 

2|-y], surface wettability p-l?!, fluid structure [7H9|, and 



interface include surface roughness 
shear rate 
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12l |. The slip length in the micron range is reported for complex flows near 
superhydrophobic surfaces 13| and flows of high molecular weight polymers [ij] , while the 
magnitude of the slip length in the range of a few tens of nanometers is typically measured 
for flows of water over smooth nonwetting surfaces [l5| . 

A number of molecular dynamics (MD) simulation studies have been carried out to in- 



vestigate the influence of structural properties o 
and flat crystalline walls on the degree of slip 16 



■ the interface between monatomic fluids 



271]. It was shown that the slip length is 



directly related to the intensity of structure induced in the first fluid layer by the periodic 
potential from the solid substrate |18i]. The slip is reduced with increasing wall- fluid inter- 
action energy and due to the formation of commensurate structures between solid wall and 
adjacent fluid layer. The solid walls are usually modeled as an array of atoms arranged on 
sites of a periodic lattice. Two types of walls are considered, solid and thermal, where the 
wall atoms are either fixed at the lattice sites or allowed to oscillate in the harmonic poten- 
tial. It was found that the slip length weakly depends on the value of the spring stiffness 
coefficient for sufficiently strong harmonic bonds |24J . |27| . In addition, the slope of the shear 
rate dependence of the slip length is not significantly affected by stiff springs [24] , except at 
very high shear rates 281 ]. 

At the interface between simple fluids and atomically smooth, weakly attractive surfaces, 
the slip length is constant at low shear rates and increases nonlinearly at higher shear 
rates, as originally discovered by Thompson and Troian 25| and later confirmed by several 



studies 
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27 



29 



30| . For sufficiently strong wall- fluid interactions and incommensurate 



structures of the liquid and solid phases at the interface, the slip length varies almost linearly 



with shear rate 
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It should be noted, however, that if the slip length at low shear rates 



is less than about a molecular diameter then the boundary conditions for dense monatomic 
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25| . Also, it was shown that molecular-scale surface 



fluids are rate-independent [18 

roughness reduces the magnitude of the slip length and the slope of its rate dependence [19, 



24 



30 



32|. 



It was recently demonstrated that the effective slip length for flows over anisotropic 
surfaces with two-component texture of different wettability is largest (smallest) for parallel 
(perpendicular) orientation of stripes with respect to the mean flow 33] • These conclusions 
hold when the stripe width is comparable to the molecular diameter 3^ . For the transverse 
orientation of the flow relative to the stripes, the slip is reduced because of the molecular 
scale corrugation of the composed surface potential, while for the parallel orientation, the 
fluid molecules are transported along homogeneous stripes with either no-slip or partial slip 
conditions, and, therefore, the effective slip length is enhanced 3^. More recently, it was 
observed that the slip length also depends on the crystal lattice plane in contact with the 
fluid and on the lattice orientation with respect to the flow direction 35|. In the present 
study, we will show that at sufficiently high shear rates, the slip flow is anisotropic for 
atomically flat crystalline surfaces; and, in particular, the slip length is enhanced when the 
shear flow is oriented along the crystallographic axis of the wall lattice. 

Recent studies of friction between adsorbed monolayers and smooth crystalline surfaces 
are relevant to the analysis of flow boundary conditions 49|, |50|. It was found that the slip 
time, which represents the transfer of momentum between the adsorbed monolayer and the 
substrate, is proportional to the phonon lifetime divided by the normalized peak value of the 
structure factor computed in the monolayer at the main reciprocal lattice vector 49|. Also, 
the simulation results have shown that the slip time is independent of the sliding direction 
if the slip velocity of the monolayer is much smaller than the speed of sound [45*] • In the 
linear regime between friction force and sliding velocity and in the range of fllm coverages 
from submonolayer to bilayer, the slip times were computed directly from the decay of the 
fllm velocity and from the decay of the velocity correlation function at equilibrium J50|. 

During the last two decades, several MD studies have examined slip boundary conditions 
at the interface between polymeric fluids and flat crystalline surfaces 36l-l46l] . The velocity 
profiles with stick boundary conditions were observed when a highly viscous interfacial layer 



was formed because of the strong wall-fluid iuteraetion energy UMM, high fluid density 

and pressure [36|, |37|, |45|, or chemical structure of chain molecules [47|. The variation of 
the slip length as a function of shear rate was reported for fiat polymer-solid interfaces with 



weak wall-fluid interactions 
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41|-|43, 



45 



46 



48| . In our previous studies J43|, |48| , it was 



shown that the rate dependence of the slip length acquires a local minimum at low shear 
rates followed by a rapid growth at higher shear rates. Furthermore, it was found that in 
a wide range of fluid densities, the friction coefficient at the liquid-solid interface undergoes 
a gradual transition from a constant value to the power-law decay as a function of the slip 
velocity 43|,|45|]. Remarkably, the data for the friction coefficient at different fluid densities 
and shear rates were found to collapse onto a master curve when plotted as a function of the 
product of the main peak in the structure factor and the contact density of the flrst fluid 
layer 43|, |45| . Although promising results were obtained, the simulations were performed at 
a single wall density and only for one orientation of the fee lattice with respect to the shear 



flow direction 



43 



45|. 



In this paper, molecular dynamics simulations were conducted for twenty systems that 
include monatomic and polymeric fluids conflned by flat crystalline surfaces. In agreement 
with previous studies, we report the nonlinear rate dependence of the slip length for systems 
with weak wall-fluid interaction energies and relatively dense walls. The simulation results 
indicate that the friction coefficient (the ratio of fluid viscosity and slip length) in the linear- 
response regime is a function of a single variable that is a product of the height of the 
normalized main peak in the structure factor and the contact density of the flrst fluid layer 
near the solid wall. We will show that the onset of the nonlinear regime between the wall 
shear stress and slip velocity is determined by the diffusion of fluid monomers within the 
flrst layer. 

The rest of the paper proceeds as follows. The details of molecular dynamics simulations 
and parameter values for twenty systems are described in the next section. The results 
for the rate-dependent shp length, friction coefficient, and fluid structure are presented in 
Section lllli The conclusions are given in the last section. 



II. MOLECULAR DYNAMICS SIMULATION MODEL AND PARAMETER VAL- 
UES 

The geometry of the computational domain and the steady flow profile are shown schemat- 
ically in Figured] The fluid undergoes planar shear flow between two atomically flat walls. 
The fluid phase consists of Nf = 9600 monomers. The interaction between any two fluid 
monomers is modeled via the truncated Lennard- Jones (LJ) potential 

'a\i2 /o-\6- 
r . 



VLj(r) = 4e 



(1) 



where e and a are the energy and length scales of the fluid phase and rc = 2.5 a is a cutoff 
radius. The interaction between the wall atoms and fluid monomers is also modeled by the 
LJ potential with parameters ^wf (listed in Tabled]) and a^{ = a. The wall atoms do not 
interact with each other via the LJ potential. 

Three types of fluid were considered in the present study, i.e., monomeric (or simple) 
fluid and polymer melts with the number of monomers per chain A^ = 10 and A^ = 20. In 
the case of polymers, the nearest-neighbor monomers in a chain interact through the finitely 



extensible nonlinear elastic (FENE) potential 51 1 

^™i.(r) = -|r„2ln[l-rVr„^], (2) 

with the energy and length parameters kg = 30 ecr^"^ and Tq = 1.5 a introduced by Kremer and 



Grest [52]]. Figure[2] shows a snapshot of an unentangled polymer melt with linear flexible 
chains A^ = 20 confined between solid walls. 

The heat exchange between the fluid phase and the external heat bath was regulated via 



a Langevin thermostat 53|, which was applied only to the direction of motion perpendicular 
to the plane of shear 18|] . The equations of motion for fluid monomers in all three directions 
are given as follows: 



mxj 



Ef^. (3) 



.,. dxi 



mjji + mTyi = - V" — ^ + fi , (4) 

^ dyi 



TnZi 



where the summation is performed over the fluid monomers and wall atoms within the cutoff 
radius rc = 2.5 a, r = 1.0r~^ is the friction coefficient, and fi is a random force with zero 
mean and variance (/i(0)/j(t)) = 2mkBTr6{t)6ij determined from the fluctuation-dissipation 
theorem. The thermostat temperature is T = 1.1 e/ks, where ks is the Boltzmann constant. 



The equations of motion were integrated using the fifth-order gear-predictor algorithm 54 1 
with a time step At = 0.002 r, where r = ynrw^/e is the characteristic time o: 



the L 
potential. The small time step At = 0.002 r was used in our previous studies [23|, |43|, |45| 
for similar MD setups in order to compute accurately the trajectories of fiuid molecules and 
wall atoms near interfaces. 

Each confining wall is composed of 1152 atoms arranged in two layers of the face-centered 
cubic (fee) or body-centered cubic (bcc) lattice. The wall density, lattice type, its orientation 
with respect to the shear fiow direction, and wall-fiuid interaction energy are reported in 
TablelB The wall atoms were either fixed at the lattice sites or were allowed to oscillate 
about their equilibrium lattice positions under the harmonic potential Vsp = | k r^ with the 
spring stiffness coefficient k = 1200 e / cx'^ . In the latter case, the Langevin thermostat was 
applied to the x, y and z components of the wall atom equations of motion. For example, 
the X component of the equation of motion is given by 

m^ Xi + m^,Txi = - ^ -— ^ - -— ^ + fi , (6) 

where rriw = 10 m, the friction coefficient is F = 1.0 r^^ and the sum is taken over the neigh- 
boring fiuid monomers within the cutoff radius r^ = 2.5 a. Periodic boundary conditions 
were imposed along the the x and y directions parallel to the confining walls. 

Initially, the fiuid was equilibrated at a constant normal pressure (shown in Table[T]) 
for about 5 x lO^r. Then, the channel height was fixed and the system was additionally 
equilibrated for 5 x lO^r at a constant density ensemble. The steady fiow was generated by 
moving the upper wall with a constant speed U in the x direction parallel to the immobile 
lower wall (see Fig.[l]). The lowest speed of the upper wall is f/= 0.05 a/r. Both fiuid 
velocity and density profiles were computed within horizontal bins of thickness A2; = 0.01 a 
for a time period up to 6 x lO^r. 



III. RESULTS 

A. Fluid density and velocity profiles 

The averaged fluid density and velocity profiles are presented in Fig.Ofor the upper wall 
speeds U = 0.5 a/r and U = 4.0 a/r. The density profiles exhibit a typical layered structure 
which extends for about 5 cr — 6 a away from the solid walls. The amplitude of the first peak 
in the density profile determines the contact density pc- We emphasize that the thickness 
of the averaging bins A2; = 0.01 a is small enough so that the magnitude of the peaks does 
not depend on the bin thickness and bin location relative to the walls. On the other hand, 
the shape of the density profiles will remain unchanged if thinner bins are used; however, 
the averaging is computationally more expensive. As evident from Fig.[3](b), the contact 
density is reduced at higher upper wall speeds. 

Two representative velocity profiles normalized by the upper wall speed are shown in 
Fig.[3](a). The slip velocity increases at higher upper wall speeds. The location of the liquid- 
solid interface (marked by the dashed vertical lines in Fig.[3]) is defined at the distance 0.5 cr 
away from the wall lattice planes to take into account the excluded volume due to wall 
atoms. The slip length was computed from the linear fit to the velocity profiles excluding 
regions of about 2 a from the solid walls. At all shear rates examined in the present study, 
the velocity profiles are linear across the channel and the slip length is larger than about 
3a. 

B. Rate dependence of viscosity and slip length 

The fluid viscosity was estimated from the relation between shear rate and shear stress 



which was computed using the Kirkwood formula 55[. The variation of viscosity as a func- 



tion of shear rate is presented in Fig. H] for selected systems listed in TablelB In agreement 



with previous studies with a similar setup 23N25|. the viscosity of monatomic fluids is inde- 



pendent of shear rate and equals /i = (2.2 ±0.2) era^^ when the fluid density is p = 0.81 cr~^. 
As expected, the shear viscosity of polymer melts with chains A^ = 10 and A^ = 20 is higher 
than the viscosity of simple monatomic fluids. For similar flow conditions, the transition 
from a Newtonian to a shear-thinning flow regime occurs at lower shear rates for polymers 
with longer chains A^ = 20 because of their slower intrinsic relaxation. The slope of the 



shear-thinning region —0.37 shown in Fig.lHis consistent with the results reported in earher 



studies for polymer melts A^ = 20 at different densities 43|, |45|. The errors arising from 
averaging over thermal fluctuations are greater at lower shear rates. 

The nonlinear rate dependence of the slip length is shown in Fig. [S] for polymer melts with 
chains A^ = 10 and A^ = 20. The shear flow direction is oriented along the crystallographic 
axis of the (111) plane of the fee wall lattice for systems 6 and 8 (see TableU]). In contrast, 
the fee lattice plane is rotated by 90° with respect to the flow direction for systems 5 and 
7 as indicated by open circles and blue vertical arrow in the inset of Fig.|5l At low shear 
rates 7r < 0.02, the slip length is independent of the wall lattice orientation relative to the 
shear flow direction; while at higher shear rates, the slip length is greater when the shear 
flow is parallel to the crystallographic axis of the triangular lattice. The same trend for the 
slip length is observed for monatomic fluids (not shown). These results demonstrate that 
at sufficiently high shear rates the slip flow is anisotropic even for atomically flat crystalline 
surfaces. 

The appearance of the local minimum in the rate dependence of the slip length reported 



in Fig. |5] for polymer melts with chains A^ = 20 was explained in the previous MD study J43|. 
The initial decay of the slip length at low shear rates is associated with a slight decrease 
in the melt viscosity while the friction coefficient at the liquid-solid interface {k = fi/Lg) 
remains constant (see also next section). With increasing shear rate, the friction coefficient 



decreases faster than the shear viscosity, and, therefore, the slip length grows rapidly 431.148 1. 
Since the transition to the shear-thinning regime occurs at higher shear rates for polymer 
melts with shorter chains A^ = 10, the slip length remains nearly constant at low shear rates 
and then increases rapidly at higher rates (see inset in Fig.|5]). These results agree well with 
the previous simulation results for slip flow of polymers with chain lengths A^ ^ 16 and 



lower fluid density 42 1. 



C. Friction coefficient versus slip velocity 

The relation between the slip length and shear rate can be expressed in terms of the 
friction coefficient at the liquid-solid interface and slip velocity. In steady-state shear flow, 
the shear stress in the bulk of the fllm (7/i) is equal to the wall shear stress (kVs). In addition, 
if the velocity proflle is linear across the channel, then by deflnition Vg = jLg and the friction 
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coefficient is given hj k= fi/Lg. In the previous MD study, the shp flow of polymer melts 
with chains A^ = 20 was studied in the range of fluid densities 0.86 ^ pa^ ^ 1.02, and 
velocity profiles were found to be linear at all shear rated examined J43|. Furthermore, the 
friction coefficient {k = fi/Lg) as a function of the slip velocity could be well fitted by the 
following equation: 

k/k* = [1 + (v;/v;)^]-°■3^ (7) 

where k* is the friction coefficient at small slip velocities when V, <^ V* and V* is the 

n 

characteristic slip velocity that determines the onset of the nonlinear regime |43[. In the 
later study, the simulations were performed at higher melt densities and the velocity profiles 
at low shear rates were curved near interfaces, and, as a result, the definition Vg = 'jLg could 



not be applied 45|]. Therefore, the friction coefficient was computed directly from the ratio 
of the wall shear stress and slip velocity of the ffist fluid layer. For polymer melt densities 
p = 1.04 cr~'^ and 1.06 cr~^, the data were also well described by Eq. ([7]), while at higher melt 



densities (1.08 ^ pa^ ^ 1.11) only the nonlinear regime was observed 45 1. 

In the present study, we extend the analysis of the friction coefficient at the interface 
between crystalline walls and polymeric fluids described by the parameters listed in TableHl 
FigureE] shows the friction coefficient as a function of the slip velocity normalized by the 
parameters k* and V* respectively. The data for all systems in Table[T] are well fitted by 
Eq. ([7]) over about three orders of magnitude. We also noticed the inverse correlation between 
the friction coefficient k* and the characteristic slip velocity V* (shown in Fig.[7]). Note that 
for every two systems with different orientation of the fee wall lattice, the values of k* are 
nearly the same, but the slip velocity V* is slightly smaller when the shear flow direction 
is parallel to the crystallographic axis. This is consistent with the dynamic response of the 
slip length reported in Fig.[5]for two different orientations of the fee lattice. 

We next argue that the onset of the nonlinear regime in Eq. ([7]) is determined by the in- 
trinsic relaxation time of the fluid monomers in the ffist layer near the solid wall. In Figure|8] 
we plot the mean square displacement of fluid monomers within the ffist layer for selected 
systems in Table[I]at equilibrium (i.e., when both walls are at rest). The displacement as a 
function of time was computed along the trajectory of a fluid monomer only if it remained 
in the first fluid layer during the time interval between successive measurements of the 
monomer position. For monatomic fluids, there is a linear dependence between the mean 
square displacement and time, and, consequently, the diffusion coefficient is well defined. 



This behavior agrees well with the exponential relaxation of the density-density autocorre- 
lation function evaluated at the wavevector of about In ja in the first layer of monatomic 



fluids confined by atomistic walls |2ll. |42|. In contrast, fluid monomers that belong to a poly- 
mer chain diffuse slower than monomers in simple fluids since their dynamics is bounded by 
diffusion of the center of mass of the polymer chain. The slope of the subdiffusive regime is 
shown in Fig. [8] by the straight dashed line. 

Finally, the comparison of the characteristic slip time of the first fluid layer and the 
diffusion time of fluid monomers between nearest minima of the surface potential is presented 
in Fig.ini The diffusion time was estimated from the mean square displacement of fluid 
monomers in the first layer at the distance between nearest minima of the periodic surface 
potential. The same distance divided by the slip velocity V* defines the characteristic slip 
time of the first fluid layer. In the case when the shear flow direction is parallel to the [110] 
fee lattice orientation (e.g., see inset in Fig.[5]), the slipping distance of the first layer was 
computed by projecting the vector, which connects nearest minima of the surface potential, 
onto the direction of flow. FigureO shows a strong correlation between the characteristic 
slip time of the adjacent fluid layer and the diffusion time of fluid monomers in that layer 
at equilibrium. These results indicate that the linear-response regime in Eq. ([7]) holds when 
the slip velocity is smaller than the diffusion velocity of fluid monomers in contact with flat 
crystalline walls. 

D. Friction coefficient and induced fluid structure 

The fluid structure near flat solid walls is characterized by density layering perpendicular 



to the surface and ordering of fluid monomers within the layers 56| . Examples of oscillatory 
density profiles in a polymer melt near confining walls were presented in Fig. [31 It is intu- 
itively expected that enhanced fluid density layering normal to the surface (obtained, for 
example, by increasing fluid pressure or wall-fluid interaction energy) would correspond to 
a larger friction coefficient at the liquid-solid interface. However, this correlation does not 
always hold; for example, the amplitude of fluid density oscillations near flat structureless 
walls might be large, but the friction coefficient is zero. As emphasized in the original paper 



by Thompson and Robbins 18|, the surface- induced fluid ordering within the first layer of 



monomers correlates well with the degree of slip at the liquid-solid interface. The measure 
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of the induced order in the adjacent fluid layer is the static structure factor, which is defined 
as follows: 






jkr, 



i=i 



(8) 



where k is a two-dimensional wavevector, Tj = {xj,yj) is the position vector of the j-th 
monomer, and A^^ is the number of monomers within the layer 18|. The probability of 
finding fiuid monomers is greater near the minima of the periodic surface potential, and, 
therefore, the structure factor typically contains a set of sharp peaks at the reciprocal lattice 
vectors. It is well established that the magnitude of the largest peak at the first reciprocal 
lattice vector is one of the main factors that determine the value of the sli p le ngth at the 
interface between fiat crystalline surfaces and monatomic fiuids [18|, l2l|, |23|, |2^ or polymer 



melts 



36 



42 



43 



45| 



Next, we discuss the infiuence of the wall-fiuid interaction energy, wall lattice type and 
orientation, and slip velocity on the structure factor computed in the first fiuid layer. The 
effect of the wall-fiuid interaction energy is illustrated in Fig. [10] for monatomic fiuids in 
contact with the (HI) plane of the fee wall lattice. The height of the surface-induced peaks in 
the structure factor is slightly larger at higher surface energy. The magnitude of the peak in 
the shear fiow direction is ^(8.65, 0) = 0.98 for e^{ = 0.3 e and ^(8.65, 0) = 1.06 for e^f = 0.4 e. 
Notice that the height of the circular ridge characteristic of short range ordering of fiuid 
monomers is larger than the amplitude of the induced peaks at the reciprocal lattice vectors. 
A similar trend in the height of the peaks in the structure factor was observed previously 
for monatomic fiuids confined by the fee walls with higher density py^ = 2.73 a^^ 23|. 

Figureim shows the structure factor computed in the first fiuid layer for polymer melts 
with chains A^ = 20 in contact with the (111) plane of the fee wall lattice. Due to the 
hexagonal symmetry of the lattice, the structure factor exhibits six peaks at the shortest 
reciprocal lattice vectors. Note that only two main peaks are present in the first quadrant. 
The magnitude of the peaks is the same at small slip velocities. The lattice orientation 
with respect to the shear fiow direction determines the location of the main peaks. Finally, 
the effect of slip velocity on the magnitude of the substrate-induced peaks in the structure 
factor is presented in Fig.[T2]for the polymer melt near the (001) plane of the bcc wall lattice. 
With increasing slip velocity, the height of the induced peak along the shear fiow direction 
decreases significantly, whereas the magnitude of the peak in the perpendicular direction is 
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less affected by slip. In contrast to monatomic ffuids, the amplitude of the peak due to short 
range order of monomers that belong to polymer chains is much smaller than the magnitude 
of the induced peaks at the shortest reciprocal lattice vectors (see Fig. fT^ . 

The correlation between surface-induced structure in the first fluid layer and the friction 
coefficient was investigate d p reviously for polymer melts with chains A^ = 20 confined by 
atomically flat walls 43|, |45|. The simulations were performed at fluid densities 0.86 ^ 
pa^ ^ 1.11 and the wall density pw = 1.40 o""^. It was found that the data for the friction 
coefficient at different shear rates and fluid densities collapsed onto a master curve when 
plotted as a function of a variable S(Oy[S{Gi) pc], where Gi is the first reciprocal lattice 



vector in the shear flow direction 



43 



45| . The collapse of the data holds at relatively small 
values of the friction coefficient A; < 4 era"^ and for slip lengths larger than approximately 
5 0". Although these results are promising, the simulations were limited to a single wall 
density and the [112] orientation of the (111) plane of the fee wall lattice. 

In the present study, a number of parameters that affect slippage at the liquid-solid 
interface have been examined, i.e., fluid and wall densities, polymer chain length, wall lattice 
type and orientation, wall-fluid interaction energy, thermal and solid walls (see Table[T]). We 
first consider the linear-response regime where the friction coefficient weakly depends on 
the slip velocity {k/k* > 0.8 in Fig.|6]). FigurefTSl shows the ratio Ls/p (an inverse friction 
coefficient) as a function of the variable 5'(0)/[5'(Gi) pd computed in the first fluid layer for 
twenty systems listed in TablelB In the case when the shear flow direction is parallel to the 
[110] orientation of the fee lattice [e.g., see Fig.[II](b)], the structure factor was computed 
at the shortest reciprocal lattice vector Gi aligned at an angle of 30° with respect to the x 
axis. Since the magnitude of the surface-induced peaks in the structure factor scales with 
the number of monomers in the first fluid layer, the height of the main peak ^(Gi) was 
normalized by the average number of monomers in the layer Ni = S{0). The data in Fig.[6] 
are well described by a power-law fit with the slope 1.13. These results suggest that at the 
interface between simple or polymeric fluids and flat crystalline surfaces, the ratio of the 
slip length and viscosity at low shear rates [ or the value of parameter k* in Eq. ([7]) ] can be 
estimated from equilibrium measurements of the structure factor and the contact density of 
the first fluid layer. 

In FigurefT^we report the dependence of the friction coefficient [k = p/Ls) on the struc- 
ture factor and contact density of the first fluid layer at all shear rates examined in this 

12 



study (the same data as in Fig. [6]). Note that at higher shear rates the derivative of Lg/ ix 
with respect to 5(0)/ [^(Gi) pc] for several systems hsted in Table[T] deviates significantly 
from the slope 1.13 shown for reference in Fig.[T3](a). In addition, for any two systems with 
the same p^ and e^f, the ratio Ls/p as a function of 5'(0)/[5'(Gi) pc] depends on the ori- 
entation of the fee wall lattice with respect to the shear flow direction. Although the data 
in Fig.[Tl](b) are somewhat scattered, the results show the same trend, namely, the friction 
coefficient decreases when the magnitude of the normalized peak in the structure factor is 

for 



reduced. The collapse of the data for Ls versus S{Gi)/ S{Q) was reported in Ref. [18 
monatomic fluids and crystalline walls when Lg < 3.5 cr and the boundary conditions are 
rate independent. In the present study, the slip lengths are greater than about 5 a except 
for the systems 13 and 16 where L^ f« 3 a at low shear rates. We flnally comment that 
our results were not analyzed with respect to the relation (between the friction coefficient 



and induced fluid structure and in-plane diffusion coefficient) derived in Ref. [2l| for simple 
fluids, because the slope of the mean square displacement versus time for polymer systems 
(shown in Fig. [8]) is less than one and thus the diffusion coefficient is not well defined. 

IV. CONCLUSIONS 

In this paper, we investigated the dynamic behavior of the slip length at interfaces be- 
tween polymeric or monatomic fluids and flat crystalline surfaces using molecular dynamics 
simulations. The polymer melt was modeled as a collection of bead-spring linear flexible 
chains below the entanglement length. We considered shear flow conditions at relatively 
low fluid densities (pressures) and weak wall-fluid interaction energies so that fluid velocity 
proflles are linear across the channel at all shear rates examined. It was found that the slip 
length does not depend on the wall lattice orientation with respect to the flow direction only 
at low shear rates, whereas the slip is enhanced at high shear rates when the flow direction 
is parallel to the crystallographic axis of the substrate. 

In the steady shear flow of either monatomic fluids or polymer melts, the friction coeffi- 
cient at the liquid-solid interface (computed from the ratio of fluid viscosity and slip length) 
undergoes a transition from a constant value to the power-law decay as a function of the 
slip velocity. The characteristic velocity of the transition is determined by the diffusion of 
fluid monomers over the distance between nearest minima of the substrate potential. It is 
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demonstrated that the friction coefficient at small slip velocities is a function of the mag- 
nitude of the surface- induced peak in the structure factor and the contact density of the 
ffist fluid layer. These conclusions hold for different wall and fluid densities, chain lengths, 
surface energies, lattice types and orientations, thermal or solid walls. 
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FIG. 1: (Color online) A schematic of the flow with slip boundary conditions at the lower and 
upper walls. Shear flow is induced by the upper wall moving with a constant speed U in the x 
direction. The slip velocity and slip length L^ are related via Vg = 'jLg, where 7 is the shear rate 
computed from the slope of the velocity profile. 
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FIG. 2: (Color online) A snapshot of fluid monomers (open blue circles) and wall atoms (filled gray 
circles) positions. Five polymer chains are marked by solid lines and filled black circles. The black 
arrow indicates the direction of the upper wall velocity U = 0.5 a /t. The fluid monomer density 
is /9 = 0.89o"^^ and the wall density is pu, = 1.80 o"^'^. The rest of parameters for the system 5 are 
given in TablelH 
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FIG. 3: (Color online) Averaged normalized velocity (a) and density (b) profiles across the channel 
for the upper wall speeds U= 0.5 a/r and U= 4.0 cj/r. The uniform monomer density of the 
polymer melt A^ = 20 away from the walls is p = 0.89(7^'^ (system 5 in Table|l]). The vertical axes 
indicate the location of the fee lattice planes (at z/a = — 11.30 and 15.14) in contact with the 
fluid. The dashed lines at z/a = — 10.80 and 14.64 denote reference planes for computing the slip 
length. 
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FIG. 4: (Color online) Shear rate dependence of the fluid viscosity fi (in units of ero"^"^) for the 
indicated systems (listed in TableUj). The dashed line with a slope —0.37 is shown for reference. 
Solid curves are a guide for the eye. 
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FIG. 5: (Color online) Slip length Lg/a as a function of shear rate for polymer melts with chains 
A^ = 20 and A^ = 10 (see inset). The system parameters are listed in TablelH Open circles in the 
inset represent the (111) face of the fee lattice atoms in contact with the fluid. The vertical blue 
arrow indicates the shear flow direction with respect to the [112] fee lattice orientation (systems 5 
and 7). The horizontal black arrow shows the flow direction with respect to the [110] orientation 
(systems 6 and 8). 
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FIG. 6: (Color online) Log-log plot of the friction coefficient k= fi/Lg (in units of era"^) as a 
function of the slip velocity Vg = Lgj (in units of cr/r) for systems listed in TablelH The values of 
the normalization parameters V* and k* are presented in Fig. [71 The dashed curves y = (l+x^)"'^'^^ 
are the best fit to the data. 
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FIG. 7: (Color online) The normalization parameters V* (in units of u/r) and k* (in units of 
eTa~'^) used to fit the data in Fig.[6]to Eq. ([7]). The indices in the inset denote systems listed in 
TableHl The dashed line with a slope —1.30 is shown for reference. 
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FIG. 8: (Color online) The mean square displacement of monomers in the first fluid layer at 
equilibrium (i.e., [/ = 0) as a function of time t (in units of r) for selected systems listed in TableHl 
The dashed lines are shown for reference. 
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FIG. 9: (Color online) A correlation between the characteristic slip time t* of the first fluid layer 
and the diffusion time t^ of fluid monomers between nearest minima of the surface potential. The 
system parameters are given in TablelH The dashed line y = x is shown as a reference. 
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FIG. 10: Two-dimensional structure factor S{kx,ky) computed in the first fluid layer for A'^ = l 
and U = 0.05 o"/r [systems (a) 19 and (b) 17 in Table[I]]. The wall-fluid interaction energy is (a) 
Ewf = 0.3 e and (b) e„f = 0.4 e. The shear flow direction (denoted by the horizontal arrow) is parallel 
to the [112] orientation of the (111) face of the fee wall lattice (open circles). 
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FIG. 11: Structure factor S{kx,ky) averaged in the first fluid layer for A^ = 20 polymer systems 
(a) 5 and (b) 6 (see parameters in TablelH). In each case, horizontal arrows indicate the shear flow 
direction with respect to the orientation of the (111) plane of the fee wall lattice (denoted by open 
circles). The upper wall speed is U = 0.05 cj/r in both cases. 
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FIG. 12: Structure factor S{kx,ky) computed in the first fluid layer for A'^ = 20 polymer system 
12 (see TableU]). The upper wall speed and slip velocity are (a) U = 0.05 cr/r and Vs = 0.012 cr/r 
and (b) U= 2.0 cr/r and Vs= 0.51 a/r, respectively. The horizontal arrows denote the shear flow 
direction with respect to the orientation of the (001) face of the bcc wall lattice (open circles). 
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FIG. 13: (Color online) Log-log plot of the ratio Ls/fJ. (in units of a'^/er) as a function of the variable 
5(0)/ [^(Gi) pc] computed in the first fluid layer at low shear rates. The system parameters are 
listed in Tablell The dashed hne y = 0.041 x^-^^ is the best fit to the data. 
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FIG. 14: (Color online) Log-log plot of the ratio Lg/fi (in units of a^/er) as a function of variables 
(a) S{0)/ [S{Gi) pc] and (b) S{0)/S{Gi) computed in the first fluid layer at all shear rates exam- 
ined. The system parameters are given in TablelH The black line with a slope 1.13 is shown for 
reference. 
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TABLE I: The fluid monomer density p, number of monomers per chain N, distance between the 
wall lattice planes in contact with fluid h, wall area in the xy plane, fluid pressure at equilibrium 
(i.e., U = 0), wall density pw, lattice type, Miller indices for the xy plane, lattice orientation along 
the shear flow direction {x direction), the x and y components of the first reciprocal lattice vector 
Gi{kx,ky), wall- fluid interaction energy, and the spring stiffness coefficient for thermal walls. 
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